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Abstract. The extraction of neutrino oscillation parameters requires the determination of the neu- 
trino energy from observations of the hadronic final state. The use of nuclear targets then requires 
the use of event generators to isolate the interesting elementary processes and to take experimental 
acceptances into account. In this talk I briefly summarize the history of event generators and their 
use in nuclear physics, talk briefly about the generators used in the neutrino community and then 
discuss future necessary developments. 
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INTRODUCTION 

The extraction of neutrino oscillation parameters from long-baseline experiments re- 
quires the knowledge of the neutrino energy. Since - due to their production mechanism 
- the neutrino beams are not monoenergetic, but instead contain a wide range of ener- 
gies, the neutrino energy has to be reconstructed from the final state of the reaction. This 
can be done by identifying a certain, fixed reaction type and using quasifree kinemat- 
ics, assuming that the reaction happens on a nucleon at rest. This is the method used by 
experiments such as MiniBooNE and K2K. Another possibility is to determine the full 
energy of the final state by calorimetric measurements; this method is used by both the 
MINOS and OPERA detectors. In case of nuclear targets the first method suffers from 
the fact that the actual reaction of a neutrino with a nucleon is not directly observable; 
final state particles, nucleons and mesons, reach the detector only after sizeable final 
state interactions (fsi). This poses a problem in, first, identifying a particular initial re- 
action (usually quasielastic (QE) scattering) and, second, in determining the final state 
energies of the first interaction from the asymptotic particles that have experienced fsi. 
In addition a natural, always -pre sent smearing due to Fermi-motion must be taken into 
account [1]. In the calorimetric method one has to deal with detection thresholds and, in 
general, detection efficiencies for all possible hadrons; the latter is also the case for the 
QE-identifying experiments that also are sensitive to detector acceptances. These known 
deficiencies in the data have to be corrected for by means of event generators [1,2] that 
introduce a sizeable dependence of the final observables on the quality of the special 
generator used [3]. A detailed comparison of predictions of generators both for QE scat- 
tering and for pion production performed for NUINT09 has revealed major differences 
between the various generators [4] . 



EVENT GENERATORS IN NUCLEAR PHYSICS 



Since the main complications in identifying elementary neutrino processes in nuclear 
targets come from nuclear final state interactions and since these are independent of 
the specific nature of the first, primary interaction it is worthwhile to look for related 
generator developments for use in quite different reactions in Nuclear Physics. Pre- 
cursors of present-day transport theoretical implementations were Monte-Carlo (MC) 
generators written to describe the final-state interactions in reactions of nucleons with 
nuclear targets. While the earliest suggestion of such methods goes back to Serber [5], 
the method has then been picked up by, among others. Metropolis et al [6], Bertini [7], 
and Cugnon [8]. About 20 years ago, with the advent of heavy-ion reactions, first meth- 
ods were developed to describe the dynamical evolution of a colliding nucleus-nucleus 
system [9-14] while taking into account the hadronic potentials and the equation of 
state of nuclear matter within the Boltzmann-Uehling-Uhlenbeck (BUU) theory. These 
codes thus went beyond the simple MC generators used until then that could not ac- 
count for effects of potentials on the propagation of particles. In addition, mesons were 
no longer instantaneously produced and, correspondingly, resonance degrees of freedom 
with their physical lifetimes were incorporated. The GiBUU code has its origins in this 
era; it started with applications to heavy-ion reactions (see [11, 13, 15] and further refer- 
ences therein). With the availability of ultra-relativistic heavy-ion beams at the AGS, the 
CERN SPS, RHIC, and the LHC this field has flourished ever since; for a more extended 
listing of transport codes presently used in heavy-ion physics see Sect. 1 in [16]. 

The first attempt to use transport-theoretical methods for the description of elementary 
processes on nuclei was an investigation of inclusive pion-nuclear reactions [17] in an 
MC calculation. A first calculation within the framework of a BUU theory of pion and 
7] photo-production off nuclei was performed by the Giessen group in [18, 19] with 
the same code that had been originally developed to describe heavy-ion collisions. Since 
then BUU theory has been used by that group to analyze a wide class of nuclear reactions 
involving - besides heavy ions - elementary projectiles such as hadrons or electrons, 
photons, and neutrinos [16]. For studies of in-medium properties of hadrons [20] it 
turned out to be essential to incorporate the transport of off-shell particles. 

All the codes and methods mentioned were developed from the outset to describe nu- 
clear phenomena. Not any special interactions were at the center of interest, but instead 
the treatment of nuclear many-body reactions and phenomena in a wide range of differ- 
ent reactions. Very independently of these developments specialized event generators for 
neutrino physics have evolved. Widely used codes such as NUANCE, NEUT, GENIE or 
new developments such as NUWRO (for a short review of neutrino generators see [21]) 
have been developed with very little coupling to the nuclear physics or the transport- 
theoretical communities. This is quite astonishing since the nuclear effects present the 
greatest challenge to the interpretation of neutrino-nucleus experiments. Only recently 
detailed studies of nuclear final state interactions, so far concentrated on n —A reactions, 
have begun [22, 23]. 



FUTURE OF GENERATORS 



In the following subsection I take up a few points of present interest that have to be 
improved in all available generators for neutrino-induced reactions on nuclei. 

Many-particle interactions. The observation by MiniBooNE that the cross section 
identified as being that of quasielastic scattering is significantly larger than the QE 
cross section obtained by various different theories and generators when they all use 
the same, world-average axial mass Ma ~ IGeV (see Fig. 15 in [24]) has triggered 
many discussions and even a call for a new paradigm [25]. Possibly connected with 
this surplus in QE scattering is the very large pion production cross section obtained in 
that experiment [26] . 

Since processes connected with initial interactions with two nucleons play a role 
in reactions with electrons on nuclei [27] it is natural to expect their influence also 
in neutrino-induced reactions. Indeed, there is now the widespread belief that the ~ 
25% surplus of events labeled as QE in the MiniBooNE experiment over theoretical 
calculations for true QE scattering is due to such so-called 2p — 2h events [28-30]. 
Thus, the widely discussed result of a significantly higher axial mass of Ma ~ 1 .2 — 1 .3 
GeV extracted from the experiments [24, 31, 32] or even Ma = 1.6 GeV in [25] could 
then simply be a consequence of trying to force a description of 2p — 2h contributions 
hy a \p — Ih model. So far none of the event generators used by the experiments to 
convert QE-like into QE events takes these 2p — 2h processes into account. It is also a 
theoretical challenge to continue the theoretical treatment of 2p — 2h contributions to 
higher energies with an aim to understand why they do not seem to contribute in the 
energy range of the NOMAD experiment [33]. 

Knock-out nucleon numbers and multi-particle interactions. Data on knock-out nu- 
cleons are directly measurable and do not require the use of a generator for the definition 
of the reaction process. Naively one expects the number of knocked-out nucleons to be 
an indication for the relative importance of QE scattering vs. other processes, since QE 
scattering (in vacuo) is connected with only one outgoing nucleon. However, due to 
fsi in a nuclear target more than one nucleon can be present in the final state even if 
the initial interaction was with one nucleon only. Fig. 1, taken from [34] and based on 
calculations with GiBUU, shows this quite well. The figure on the left shows the cross 
sections for multiparticle knock-out if the initial interaction is only QE scattering (which 
is not directly observable on a nuclear target). The figure on the right shows the actually 
observable cross section for 1 . . .4 nucleons. The cross section for seeing more than one 
nucleon in the final state increases as a function of neutrino energy and the probability 
to see 2 nucleons amounts to 20% of the 1 -nucleon contribution at 1 GeV even though 
the initial, first interaction always happened on one nucleon only. 

It is to be expected that this number will increase if the 2p — 2h contributions that can 
amount to about 25% in the total cross section are taken into account in the generators. It 
is clearly necessary to implement such processes explicitly in any event generator since 
the many particle final states will affect the energy reconstruction. The recent proposal 
by Bodek et al. [35] to include such processes simply by increasing the single-particle 
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FIGURE 1. (Color online) Cross section for knock-out of various numbers of nucleons from a ^^C 
target as a function of neutrino energy, calculated with GiBUU. The picture on the left shows the result 
obtained from primary QE events only, the right picture contains the result obtained from all primary 
interactions (from [34]). 



transverse response misses this important physics since it leads to incorrect knock-out 
nucleon numbers in the final state of the initial interaction. 

The latter is also true when treating effects of groundstate correlations [36, 37] in the 
impulse approximation (lA) where the first interaction happens on one (dressed) nucleon 
only. The spectral functions used there contain the effects of these correlations in terms 
of nontrivial connections between energy and momentum of the nucleons, i.e. spectral 
functions. There are, however, many-body contributions, for example those connected 
with exchange currents, that cannot be taken care off by an lA using single particle 
spectral functions. These have to be added explicitly into event generators. 

For handling these off-shell effects in the spectral functions in the very first, primary 
interaction there are methods (recipes) available that have been developed over the last 
30 years [38, 39]; these methods seem to be accepted in the community, but suffer 
from problems connected with insufficient knowledge of off-shell properties of the 
interaction vertices and gauge invariance. By assumption they furthermore do not allow 
for production and further propagation of off-shell nucleons. This restriction is hard 
to justify in the presence of fsi. Based on the seminal work of Kadanoff and Baym 
[40] developed around 40 years ago practical implementations of transport of off-shell 
particles have become available only during the last decade. While GiBUU does contain 
this new physics, all other transport codes used in neutrino physics do not. That off- 
shell transport can make a difference for knock-out nucleon cross sections and particle 
production close to threshold was shown in [41, 42]. 

Broad-Band Generators. Various reaction mechanisms, such as true QE, pion pro- 
duction, resonance excitation, DIS and 2p — 2h contributions dominate in different 
regimes of the neutrino energy so that any broad-band experiment necessarily contains 
contributions from different reaction mechanisms. To isolate quasielastic scattering, for 
example, for use in the neutrino energy reconstruction then requires the use of event 
generators and theories which are reliable not only for QE scattering, but are well tested 
and reliable for a broad class of relevant reactions. In particular, the errors connected 



with the use of these generators have to be well under control [2, 3]. The latter can be 
assessed only by comparison with results from different experiments. 

Verification of Generators. There are lots of data about photo- and electro- 
production of mesons on nuclear targets, in the energy range of a few 100 MeV 
up to 200 GeV which cover the energy regime relevant for present day's neutrino exper- 
iments. These reactions involve a closely related primary interaction and the very same 
fsi as in the neutrino-induced reactions. They could thus be used to check the reliability 
and accuracy of the generators used. While it seems that some of such checks were 
done (but not published) with NEUT [43] none have been done with the other standard 
neutrino generators; on the contrary, GiBUU was developed to handle such interactions 
[16]. The necessary code changes are minimal (only the initial axial coupling has to be 
turned off). Such comparisons are ultimately more meaningful than just comparisons of 
generators with other generators. 

SUMMARY 

Mainly because of the presence of high-statistics data obtained during the last few 
years at K2K, MiniBooNE (and soon T2K, Minerva and NOvA) our understanding of 
reactions of neutrinos with nuclei has improved a lot. It seems to be quite clear now 
that the initial analyses that led to very different, higher values of the axial mass than the 
accepted world-average value of around 1 GeV were incorrect because they tried to force 
a description of 2p — 2h effects in the primary interaction in terms of a one-body model. 
Obvious is now that all generators have to be extended to include the 2p — 2h effects; 
what is open is still the question if these effects are unambiguously observable and how 
their inclusion will affect the energy reconstruction. Since these two-body couplings 
are well known from interactions of real or virtual photons with nuclei the lesson to 
be learned from this recent development is that neutrino generators do not only have 
to be compared with each other. Foremost, they have to be compared with electron- or 
photon-induced data on nuclear targets since here the primary interaction is quite similar 
to that for neutrinos and the final state interactions are the same. Very few, if none, of 
such comparisons exist for the 'working-horse generators'. 

Finally, it is clear that the extraction of neutrino oscillation parameters depends to a 
large extent on the reliability and quality of the event generators used. The quality of 
these generators could only benefit from a much closer contact to nuclear theory and 
to other generator developments in nuclear physics. This is still a long way to go for 
neutrino generators until they reach the level of sophistication that the experimental 
equipment in these long baseline experiments already has. 
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